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Abstract. In the framework of the theory of open systems based on completely 
positive quantum dynamical semigroups, we give a description of the dynamics of 
. entanglement for a system consisting of two uncoupled harmonic oscillators interacting 

with a thermal environment. Using Peres-Simon necessary and sufficient criterion for 
separability of two-mode Gaussian states, we describe the evolution of entanglement 
, in terms of the covariance matrix for a Gaussian input state. For some values 

^ , ' of the temperature of environment, the state keeps for all times its initial type: 

separable or entangled. In other cases, entanglement generation, entanglement sudden 
death or a repeated collapse and revival of entanglement take place. We determine 
the asymptotic Gaussian maximally entangled mixed states (GMEMS) and their 
. corresponding asymptotic maximal logarithmic negativity. 
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1. Introduction 

Quantum entanglement of Gaussian states constitutes a fundamental resource in 
continuous variable quantum information processing and communication [1]. For the 
class of Gaussian states there exist necessary and sufficient criteria of entanglement [21 [3] 
and quantitative entanglement measures [U [5]. The quantum information processing 
tasks are difficult to implement, due to the fact that any realistic quantum system 
is not isolated and it always interacts with its environment. Quantum coherence and 
entanglement of quantum systems are inevitably influenced during their interaction with 
the external environment. As a result of the irreversible and uncontrollable phenomenon 
of quantum decoherence, the purity and entanglement of quantum states are in most 
cases degraded. Therefore in order to describe realistically quantum information 
processes it is necessary to take decoherence and dissipation into consideration. 
Decoherence and dynamics of quantum entanglement in continuous variable open 
systems have been intensively studied in the last years [HI [TJ El [HJ [THl ttH [ISl [13] . 

When two systems are immersed in an environment, then, in addition to and 
at the same time with the quantum decoherence phenomenon, the environment can 
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also generate a quantum entanglement of the two systems and therefore an additional 
mechanism to correlate them |9l [HI [15] . In this paper we study, in the framework of the 
theory of open systems based on completely positive quantum dynamical semigroups, the 
dynamics of the continuous variable entanglement of two identical harmonic oscillators 
coupled to a common thermal environment. We are interested in discussing the 
correlation effect of the environment, therefore we assume that the two oscillators are 
uncoupled, i.e. they do not interact directly. The initial state of the subsystem is 
taken of Gaussian form and the evolution under the quantum dynamical semigroup 
assures the preservation in time of the Gaussian form of the state. In section 2 
we write the Markovian master equation in the Heisenberg representation for two 
uncoupled harmonic oscillators interacting with a general environment and the evolution 
equation for the covariance matrix of the considered subsystem. By using the Peres- 
Simon criterion for separabihty of two- mode Gaussian states [21 |16], we investigate 
in section 3 the dynamics of entanglement for this system. For certain values of the 
environment temperature, the state keeps for all times its initial type: separable or 
entangled. For other values of the temperature, entanglement generation, entanglement 
sudden death or a repeated collapse and revival of entanglement take place. In recent 
papers, in a different model, Paz and Roncaglia [T71 [18] studied numerically the exact 
entanglement behaviour of two identical oscillators in an infinite bath of other oscillators 
by using the exact master equation for quantum Brownian motion and showed that 
the entanglement can undergo three qualitatively different dynamical phases: sudden 
death, sudden death and revival, and no sudden death of entanglement. Finally, we 
determine the asymptotic Gaussian maximally entangled mixed states (GMEMS) and 
the corresponding asymptotic maximal logarithmic negativity, which characterizes the 
degree of entanglement of these states. A summary is given in section 4. 

2. Time evolution of covariance matrix for two harmonic oscillators 

We study the dynamics of the subsystem composed of two identical non-interacting 
oscillators in weak interaction with a thermal environment. In the axiomatic formalism 
based on completely positive quantum dynamical semigroups, the irreversible time 
evolution of an open system is described by the following general quantum Markovian 
master equation for an operator A in the Heisenberg representation (f denotes Hermitian 
conjugation) [HI [20]: 

f = ^[H, A] + ^ E(2^'^^. - ViV.A - AV^V,). (1) 

j 

Here, H denotes the Hamiltonian of the open system and the operators Vj, Vj , defined 
on the Hilbert space of H, represent the interaction of the open system with the 
environment. 

We are interested in the set of Gaussian states, therefore we introduce such quantum 
dynamical semigroups that preserve this set during time evolution of the system. 
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Consequently H is taken to be a polynomial of second degree in the coordinates x, y 
and momenta Px,Py of the oscillators and Vj,Vj are taken polynomials of first degree 
in these canonical observables. Then in the linear space spanned by coordinates and 

ED: 



momenta there exist only four linearly independent operators Vj=i.2,3,4 



(2) 



Xj , dyj , Uxj , Uyj 



bxj, byj are complex coefficients. The Hamiltonian H of the two uncoupled 



where a. 

identical harmonic oscillators of mass m and frequency u is given by 



H 



-iPx+Pt) + ^ix' + y') 



(3) 
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The fact that the evolution is given by a dynamical semigroup implies the positivity 
of the following matrix formed by the scalar products of the four vectors ax,ay,bx, by 

respectively: 

/ (axa^) (a^b^) (a^ay) (a^by) \ 



whose entries are the components a^j, ayj, b^j, byj, 



(bxax) 



[cLxbx) 
(bxbx) 
{aybx) 
(bybx) 



[axO-y) 
(b^ay) 
{aytty) 
(byay) 



(bxby) 



{ttyby] 

V (bya^) (byK) (byay) {byby j 

We take this matrix of the following form, where all diffusion coefficients D 
and dissipation constant A are real quantities (we put from now on ^ = 1 



(4) 



XXI Dxp^ )• • ■ 



/ 


Dxx 


—D — i- 

^xpx i 2 


Dxy 


Dxpy 








^PxPx 


~Dyp^ 


DpxPy 






Dxy 


~^ypx 


Dyy 


~^ypv ~ 


i^ 


V 


~Dxpy 


^PxPy 


~Dypy + ig 


DpyPy 


/ 



(5) 



It follows that the principal minors of this matrix are positive or zero. From the Cauchy- 
Schwarz inequality the following relations hold for the coefficients defined in ([5]): 



DxxDp^p^ 



D 



XPx 



)2 

— ^yy^PyPv 



D^^Dyy - D% > 0, 



DxxDp^p^ - Dt^^ > 0, 



DpxPx DpyPy 

DyyDp^p^ 



X2 

>- 

ypy — 4 ' 

Kpy > 0' 

- dL > 0. 



■ypx - - (6) 
A two-mode Gaussian state is completely characterized by its first and second 
moments of canonical variables. We therefore introduce the following 4x4 bimodal 
covariance matrix, which entirely specifies a two-mode Gaussian state (all first moments 
have been set to zero by means of local unitary operations which do not affect the 
entanglement): 



ait) 



( (7xx{t) 



Or 



a. 



<^r,^Vx{t) 



cr. 



xy 



(^xy{t) 
\(^xpy{t) 



'PxPx 

^ypx 



0"„r,^ (t) 



it) 

(^ypxit) 
a„y{t) 



cr. 



PxPy 



'yy\ 
'ypy^ 



(yim.i{t) 



xPy 
PxPy 

ypy 

PyPy 



(Jn^nyif) 



a. 



it) 



(7) 



crr,.,„.,(t) / 



The problem of solving the master equation for the operators in Heisenberg 
representation can be transformed into a problem of solving first-order in time, coupled 
linear differential equations for the covariance matrix elements. Namely, from ([T]) we 
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obtain the following system of equations for the quantum correlations of the canonical 
observables, written in matrix form |2T] (T denotes a transposed matrix): 

da{t) 



dt 



Ya{t) + a{t)Y' +2D, 



where 



/ -A 



Y 



-mu 





1/m 
-A 












-A 



—muj 



\ 


1/m 
-A / 



(8) 



(9) 



D 



D 



D 
D 



PxPx 



D 
D 



xy 



ypx 



Dxpy 
D. 



\ 



PxPy 



D 



xy 



\D 



XPy 



D 
D 



ypx 



PxPy 



D 
D 



yy 



D 



ypy 



(10) 



ypy 



^PyPy / 



The time-dependent solution of Eq. ([8]) is given by [21] 

a{t) = M{t)[a{0) - a{oo)]M^{t) + a{oo), (11) 

where the matrix M{t) = exp(yt) has to fulfill the condition limt^oo M{t) = 0. In order 
that this limit exists, Y must only have eigenvalues with negative real parts. The values 
at infinity are obtained from the equation 

Ya{oo) + a{oo)Y^ = -2D. (12) 



3. Dynamics of two-mode continuous variable entanglement 

The characterization of the separability of continuous variable states using second-order 
moments of quadrature operators was given in Refs. [2, 3J. A two-mode Gaussian 
state is separable if and only if the partial transpose of its density matrix is non- 
negative [necessary and sufficient positive partial transpose (PPT) criterion]. A two- 
mode Gaussian state is entirely specified by its covariance matrix ([7]), which is a real, 
symmetric and positive matrix with the block structure 

where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric 
covariance matrices for the individual one-mode states, while the matrix C contains 
the cross-correlations between modes. Simon |2] derived a PPT criterion for 
bipartite Gaussian continuous variable states: the necessary and sufficient criterion for 
separability is S{t) > 0, where 

Sit) = det A det 5 + Q - I det 

-Tt[AJCJBJC^J] - ^ (det A + det B) (14) 
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and J is the 2x2 symplectic matrix 

Since the two oscillators are identical, it is natural to consider environments for which 
Dxx = Dyy, Dxp^ = Dyp^, Dp^p^ = D p^p^ , D ^Py = Dyp^. TYiQii bothunimodal covariance 
matrices are equal, A = B, and the entanglement matrix C is symmetric. 

3.1. Time evolution of entanglement 

In order to describe the dynamics of entanglement, we use the PPT criterion [21 [16] 
according to which a state is entangled if and only if the operation of partial transposition 
does not preserve its positivity. Concretely, we have to analyze the time evolution of the 
Simon function S{t) (fT^ . For a thermal environment characterized by the temperature 
T, we consider such diffusion coefficients, for which 

muDxx = ^ = ^coth^, D,p^ = 0, (16) 



moo 



m u Dxy = Dp^p^. (17) 

This corresponds to the case when the asymptotic state is a Gibbs state [20]. We 
consider two cases, according to the type of the initial Gaussian state: 1) separable and 
2) entangled. 

1) To illustrate a possible generation of the entanglement, we represent in Figure 
1 the dependence of function S{t) on time t and temperature T for a separable initial 
Gaussian mixed state. We notice that, according to Peres-Simon criterion, for relatively 
small values of the temperature T, the initial separable state {S{t) = 0) becomes 
entangled shortly after the initial moment of time t = 0. For relatively large values 
of T, S{t) is strict positive and the state remains separable for all times. 

Depending on the environment temperature, there are three situations in the case 
of a generated entanglement [22]: a) entanglement may persist forever, including the 
asymptotic final state; b) there exist repeated collapse and revival of entanglement; 
c) the entanglement is created only for a short time, then it disappears and the state 
becomes again separable. The entanglement of the two modes can be generated from 
an initial separable state during the interaction with the environment only for certain 
values of mixed diffusion coefficient D^py and dissipation constant A. 

2) An example of the evolution of an entangled initial state is illustrated in Figure 
2, where we represent the dependence of function S{t) on time t and temperature T for 
an entangled initial Gaussian mixed state. For relatively small values of T, the initial 
entangled state may remain entangled for all times [22]. For relatively large values 
of temperature T, at some finite moment of time, S{t) takes non-negative values and 
therefore the state becomes separable. This is the so-called phenomenon of entanglement 
sudden death. This phenomenon is in contrast to the loss of quantum coherence, which 
is usually gradual [12]. Depending on the values of the temperature, it is also possible 
to have a repeated collapse and revival of the entanglement. 
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Figure 1. Dependence of Simon separability function S on time t and environment 
temperature C = coth for A = 0.1, D^p^ = 0.049 and for a separable initial 
Gaussian mixed state with initial correlations axxiO) — I7 cTp^pxCO) = 1/2, f^xp^i^^) = 
cTxyiO) = <^p^py{0) = o-xp„(0) = 0. We take m = uj = h^l. 




Figure 2. Same as in Figure 1, for an entangled initial Gaussian mixed state with 
initial correlations cr^xiO) = 1, o-p=„px(0) = 1/2, crxp,(0) = 0, a^yiO) = 1/2, crp,p„(0) = 
-1/2, a,p„(0) =0. 



3.2. Logarithmic negativity 

For Gaussian states, the measures of entanglement of bipartite systems are based 
on some invariants constructed from the elements of the covariance matrix [6l ITU] . 
Heisenberg uncertainty principle can be expressed as a constraint on the global 
symplectic invariants A = det A + det B + 2 det C and det a J2]: 

A<7 + 4detcr, (18) 
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and the symplectic eigenvalues u^, which form the symplectic spectrum of covariance 
matrix a are determined by 



2z/2 = A T VA2 - 4 det a. (19) 

In terms of relation ( IT8|) takes the form z/_ > 1/2. 

In order to quantify the degrees of entanglement of the infinite-dimensional bipartite 
system states of the two oscillators it is suitable to use the logarithmic negativity. For 
a Gaussian density operator, the logarithmic negativity is completely defined by the 
symplectic spectrum of the partial transpose of the covariance matrix. It is given by 
En = max{0, — log2 2z/_}, where i>_ is the smallest of the two symplectic eigenvalues of 
the partial transpose a of the 2-mode covariance matrix a : 



2z>2 = At V^^ -4det(T. (20) 

Here A is given by A = A — 4 det C = det A + det B — 2 det C. A state is separable if 
and only if 

i>_ > 1/2 (21) 

and logarithmic negativity quantifies the violation of inequality (I2ip . 
In our model, the logarithmic negativity is calculated as 

E^(t) = -ilog2[4/(a(t))], (22) 



where 



f{a{t)) = ^(det A + det B) - det C 

■ I -,2 \ 1/2 



2 (det A + det B) - det C 



det a{t) . (23) 



It determines the strength of entanglement for E^^t) > 0, and if E^it) < 0, then the 
state is separable. 

As expected, the logarithmic negativity has a behaviour similar to that one of the 
Simon function in what concerns the characteristics of the state of being separable or 
entangled [23l El [251 [26] . 

3.3. Asymptotic entanglement 

From (fT2l) . ( [T6l) and (fTTj) we obtain the following elements of the asymptotic matrices 
A{oo) = B{oo) : 

mua^^ioo) = ^E^P^ = 1 coth a,p^{oo) = (24) 

mu 2 2ki 

and of the entanglement matrix C(oo) : 

°« = ^ Ji^2l. 2^ 25 

m^A(A^ T w ) 



(^xp,(oo) = ayp,(oo) = 4-?^' (26) 
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The mixedness of a quantum state p is characterized by its purity /i = Trp^. 
For a two- mode Gaussian state the purity is given by yU = l/Ay/deta. The marginal 
purities /ij (z = 1, 2) of the reduced states in mode i are given by ni = 1/2V det A and 
H2 = 1/2 V det B. The global and marginal purities range from to 1, and they fulfill 
the constraint /i > /ii/i2, as a direct consequence of Heisenberg uncertainty relations. In 
Ref. [27] the following upper and lower bounds on the invariant A have been obtained 
in terms of global and marginal purities: 

2/i 4/if/i^ 

The invariant A has a direct physical interpretation |27]: at given global and marginal 
purities, A determines the amount of entanglement of the state and the smallest 
symplectic eigenvalue U- of the partially transposed state is strictly monotone in A. 
Consequently, the entanglement of a Gaussian state with fixed global purity /i and 
marginal purities pi, fi2 is strictly increasing with decreasing A. According to double 
inequality f l25]) . giving lower and upper bounds on A, there exist both maximally and 
minimally entangled Gaussian states. 

We analyze the existence of the entanglement in the asymptotic regime in the 
symmetric situation A = B. First we consider the particular case D^y = 0. In the limit 
of long times, we obtain the following quantities: 

detA = detB = ^, detC = - — , det a = - — , (29) 



4 ' A2 



Ct = coth d = D^p. = uj^ + A^. (30) 



where we used the notations 

u 

Then we obtain 

l = C|-4^, — = — = Ct, A = 2Vd^=^. (31) 
/i A"' /ii /i2 2/i 

These values saturate the lower bound in inequalities (l28l) and this situation entails a 
maximal entanglement. Consequently, the corresponding states are Gaussian maximally 
entangled mixed states (GMEMS). They are thermal squeezed states with the squeezing 
parameter r given by tanh2r = (I/ACt- In the pure case (det a = 1/16), these states 
are equivalent to two-mode squeezed vacua with the squeezing parameter determined 



only by the temperature of the thermal environment : tanh2r = — 1/2Ct. 
According to Ref. [27], these states are separable in the range 

/X < . (32) 
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Then, for a given temperature T, we obtain that the asymptotic final state is separable 
for the following range of positive values of the mixed diffusion coefficient d: 
Id 

^ < Ct - 1. (33) 

We remind that, according to inequalities (E]), the coefficients have to fulfill also the 
constraint ACt/2 > d. For a given temperature of the environment and for this range 
of mixed diffusion coefficients, no entanglement can occur for these states. Outside this 
region, i.e. for a temperature and diffusion coefficient satisfying 

(^T - 1 < ^ < Cr + 1, (34) 

they are Gaussian maximally entangled mixed states (GMEMS). 
The asymptotic logarithmic negativity has the form 



E^(oo) = -log2 Ct-— . (35) 




Outside the separable region, this is the maximum possible value of the logarithmic 
negativity, attained by GMEMS. It depends only on the mixed diffusion coefficient, 
dissipation constant and temperature, and does not depend on the initial Gaussian 
state. 

When (i = 0, but D^y 7^ 0, then detC > and the asymptotic final state is 
always separable. If both these diffusion coefficients are non-zero, ci 7^ 0, D.j,y ^ 0, then 
the lower bound in inequalities f l28|) is not anymore saturated and, consequently, the 
corresponding states are Gaussian non-maximally entangled mixed states. 

We notice that the asymptotic persistence of entanglement is the result of the 
competition between thermal decoherence, determined by the uni-modal diffusion 
coefficients f|T6l) and the statistical coupling of the two modes, due to the presence 
of position-position and momentum-momentum diffusion coefficients, which are chosen 
of equal relevance in our model {w?uj'^Dxy = Dp^p^) and position- momentum diffusion 
coefficients Dr^p^ = Dyp^ = d. 



4. Summary 

In the framework of the theory of open quantum systems based on completely 
positive quantum dynamical semigroups, we investigated the Markovian dynamics of 
the quantum entanglement for a subsystem composed of two noninteracting modes 
embedded in a thermal environment. By using the Peres-Simon necessary and sufficient 
criterion for separability of two-mode Gaussian states, we have described the evolution of 
entanglement in terms of the covariance matrix for Gaussian input states. The dynamics 
of the quantum entanglement is sensitive to the initial states and the parameters 
characterizing the environment (diffusion and dissipation coefficients and temperature). 
For some values of these parameters, the state keeps for all times its initial type: 
separable or entangled. In other cases, entanglement generation or entanglement 
suppression (entanglement sudden death) take place or even one can notice repeated 
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collapse and revival of entanglement. We have also shown that, independent of the 
type of the initial state - separable or entangled, for certain values of temperature, the 
initial state evolves asymptotically to an equilibrium state which is entangled, while for 
other values of temperature the asymptotic state is separable. For a given temperature, 
we calculated the range of mixed diffusion coefficients which determine the existence of 
asymptotic Gaussian maximally entangled mixed states (GMEMS). We obtained also 
the expression of the maximal logarithmic negativity, which characterizes the degree of 
entanglement of these states. 
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